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Abstract We review some recent results of the fractional variational calculus. Nec- 
essary optimality conditions of Euler-Lagrange type for functionals with a La- 
grangian containing left and right Caputo derivatives are given. Several problems 
are considered: with fixed or free boundary conditions, and in presence of integral 
constraints that also depend on Caputo derivatives. 
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1 Introduction 

Fractional calculus plays an important role in many different areas, and has proven 
to be a truly multidisciplinary subject [20,30]. It is a mathematical field as old as the 
calculus itself. In a letter dated 30th September 1695, Leibniz posed the following 
question to LHopital: "Can the meaning of derivative be generalized to derivatives 
of non-integer order?" Since then, several mathematicians had investigated Leib- 
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niz's challenge, prominent among them were Liouville, Riemann, Weyl, and Let- 
nikov. There are many applications of fractional calculus, e.g., in viscoelasticity, 
electrochemistry, diffusion processes, control theory, heat conduction, electricity, 
mechanics, chaos and fractals, and signals and systems [12,22]. 

Several methods to solve fractional differential equations are available, using 
Laplace and Fourier transforms, truncated Taylor series, and numerical approxima- 
tions. In [7] a new direct method to find exact solutions of fractional variational 
problems is proposed, based on a simple but powerful idea introduced by Leitmann, 
that does not involve solving (fractional) differential equations [32]. By an appro- 
priate coordinate transformation, we rewrite the initial problem to an equivalent 
simpler one; knowing the solution for the new equivalent problem, and since there 
exists an one-to-one correspondence between the minimizers (or maximizers) of the 
new problem with the ones of the original, we determine the desired solution. For a 
modern account on Leitmann's direct method see [25,26]. 

The calculus of variations is a field of mathematics that deals with extremiz- 
ing functionals [33]. The variational functionals are often formed as definite inte- 
grals involving unknown functions and their derivatives. The fundamental problem 
consists to find functions y(x), x £ [a,b], that extremize a given functional when 
subject to boundary conditions y(a) = y a and y(b) = yt,. Since this can be a hard 
task, one wishes to study necessary and sufficient optimality conditions. The sim- 
plest example is the following one: what is the shape of the curve y(x), x £ [a,b], 
joining two fixed points y a and yi,, that has the minimum possible length? The an- 
swer is obviously the straight line joining y a and yt,. One can obtain it solving the 
corresponding Euler-Lagrange necessary optimality condition. If the boundary con- 
dition y(b) = yt, is not fixed, i.e., if we are only interested in the minimum length, 
the answer is the horizontal straight line y(x) = y a , x £ [a,b] (free endpoint prob- 
lem). In this case we need to complement the Euler-Lagrange equation with an 
appropriate natural boundary condition. For a general account on Euler-Lagrange 
equations and natural boundary conditions, we refer the reader to [23,24] and refer- 
ences therein. Another important family of variational problems is the isoperimetric 
one [5]. The classical isoperimetric problem consists to find a continuously differ- 
entiable function y = y(x), x £ [a,b], satisfying given boundary conditions y(a) = y a 
and y(b) = yt,, which minimizes (or maximizes) a functional 



The most famous isoperimetric problem can be posed as follows. Amongst all closed 
curves with a given length, which one encloses the largest area? The answer, as we 
know, is the circle. The general method to solve such problems involves an Euler- 
Lagrange equation obtained via the concept of Lagrange multiplier (see, e.g., [4]). 




subject to the constraint 
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The fractional calculus of variations is a recent field, initiated in 1997, where 
classical variational problems are considered but in presence of some fractional 
derivative or fractional integral [31]. In the past few years an increasing of interest 
has been put on finding necessary conditions of optimality for variational problems 
with Lagrangians involving fractional derivatives [1,9-1 1, 16-19,27, 28], fractional 
derivatives and fractional integrals [3,6, 15], classical and fractional derivatives [29], 
as well as fractional difference operators [13, 14]. A good introduction to the sub- 
ject is given in the monograph [21]. Here we consider unconstrained and constrained 
fractional variational problems via Caputo operators. 



2 Preliminaries and notations 

There exist several definitions of fractional derivatives and fractional integrals, 
e.g., Riemann-Liouville, Caputo, Riesz, Riesz-Caputo, Weyl, Grunwald-Letnikov, 
Hadamard, and Chen. Here we review only some basic features of Caputo's frac- 
tional derivative. For proofs and more on the subject, we refer the reader to [20, 30]. 

Let / : [a, b] — > M. be an integrable function, a > 0, and T be the Euler gamma 
function. The left and right Riemann-Liouville fractional integral operators of order 
a are defined by 2 



and 



'(«)■ 

respectively. The left and right Riemann-Liouville fractional derivative operators of 
order a are, respectively, defined by 



d n 



and 

d" 

n a — C-1 V 7 -— - o l"~ a 
x b - ( lj dx" Jh ' 

where n = [a] + 1 . Interchanging the composition of operators in the definition of 
Riemann-Liouville fractional derivatives, we obtain the left and right Caputo frac- 
tional derivatives of order a: 

dx" 



2 Along the work we use round brackets for the arguments of functions, and square brackets for the 
arguments of operators. By definition, an operator receives a function and returns another function. 
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and 

c x Dl := x I" b - a ° (-!)"£• 



Theorem 1. Assume that f is of class C n on [a,b]. Then its left and right Caputo 
derivatives are continuous on the closed interval [a,b]. 

One of the most important results for the proof of necessary optimality con- 
ditions, is the integration by parts formula. For Caputo derivatives the following 
relations hold. 

Theorem 2. Let a > 0, and /,g : [a,£>] — > R foe C" functions. Then, 
f b g(x) ■ c a D a x [/] (x)dx = ["fix) ■ x D a b [g] {x)dx 

+ tlDt +] -'\g]{x)- x D'l- X - j [f]{x)\ b 
j=o l ja 

and 

t g(x) ■ c x D a b [/] (x)dx = ("fix) ■ a D« [g] (x)dx 

J a J a 

7=0 

where a D x = a I~ k and x D b = X I7 whenever k < 0. 

In the particular case when < a < 1, we get from Theorem 2 that 

( b g{x)- C a D"[f]{x)dx = tf{x) ■ x D«[g](x)dx + [4- a [g](x) .f(x)] b a 

J a J a 

and 

fg(x) - c x D a b [f](x)dx = f f(x)- a D«[g](x)dx- [J x - a [g](x) ■f(x)] b a . 

J a J a 

In addition, if / is such that f{a) = f(b) = 0, then 

g{x)- c a D a x [f}(x)dx = f 'f{x) -xDtlg^dx 

J a 



and ^ ^ 

' g(x)- c x D?[f](x)dx= [ f(x)- a D^[g](x)dx. 

J a 

Along the work, we denote by <3,L, i = l,...,m (m 6 N), the partial derivative of 
function L : W" — > M. with respect to its ith argument. For convenience of notation, 
we introduce the operator c a [■] a defined by 
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c a [y]p :=^(i,y(i),^MW,^MW) - 
where a,j3 G (0,1). 



3 Euler-Lagrange equations 

The fundamental problem of the fractional calculus of variations is addressed in the 
following way: find functions y G #, 

:= {.v G C 1 ( [a, 6] ) | y (a) = y a and y(b) =y b }, 

that maximize or minimize the functional 

J(y)= [ (Lo c a \y} p )(x)dx. (1) 

J a 

As usual, the Lagrange function L is assumed to be of class C 1 on all its arguments. 
We also assume that d^Lo^y]^ has continuous right Riemann-Liouville fractional 
derivative of order a and AiZo^Jy]^ has continuous left Riemann-Liouville frac- 
tional derivative of order /3 for y G <S . 

In [1] a necessary condition of optimality for such functionals is proved. We 
remark that although functional (1) contains only Caputo fractional derivatives, 
the fractional Euler-Lagrange equation also contains Riemann-Liouville fractional 
derivatives. 

Theorem 3 (Euler-Lagrange equation for (1)). Ify is a minimizer or a maximizer 
of J on S, then y is a solution of the fractional differential equation 

(d 2 Lo C a [y] p ) (X)+ X D% [d 3 Lo C a [y] p ] (x)+aPZ [d 4 Lo C a [y] p ] (x) = (2) 

for all x G [a,b]. 

Proof. Given |e| ^ 1, consider h G V where 

V := {h£C 1 ([a,b])\h(a)=0andh(b)=0}, 
and a variation of function y of type y + eh. Define the real valued function j(e) by 

rb 

j(e) = J(y + eh) = {Lo c a [ y + eh] p ) {x)dx. 

J a 

Since £ = is a minimizer or a maximizer of j, we have /(0) = 0. Thus, 
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rb 



{d 2 Lo C a {y]p) (x)-h(x) + (d 3 Lo C a \y} p ) {x)- C a D a x [h\{x) 



(d 4 Lo c a [y] p )(x)- C x D P b [h](x) 



dx = 0. 



We obtain equality (2) integrating by parts and applying the classical fundamental 
lemma of the calculus of variations [33]. 

We remark that when a — > 1, then (1) is reduced to a classical functional 

rb 



J{y)= [ f(x,y(x),y'(x))dx, 



and the fractional Euler-Lagrange equation (2) gives the standard one: 

d 2 f{x,y{x),y'(x)) - ^-d 3 f(x,y(x),y'(x)) = 0. 
dx 

Solutions to equation (2) are said to be extremals of (1). 



4 The isoperimetric problem 

The fractional isoperimetric problem is stated in the following way: find the mini- 
mizers or maximizers of functional J as in (1), over all functions y G $ satisfying 
the fractional integral constraint 

r b 



i(y)= I (8° c a \y]is)(x)dx = i. 

J a 



Similarly as L, g is assumed to be of class C 1 with respect to all its arguments, func- 
tion d^go a[y]p is assumed to have continuous right Riemann-Liouville fractional 
derivative of order a and Atgo^[y]o continuous left Riemann-Liouville fractional 
derivative of order j3 for y £ S. A necessary optimality condition for the fractional 
isoperimetric problem is given in [8]. 

Theorem 4. Let y be a minimizer or maximizer of J on S, when restricted to the set 
of functions z G S such that I(z) = L In addition, assume that y is not an extremal 
of I. Then, there exists a constant X such that y is a solution of 

(dlF o c a [y] p ) (x) + X D« [d 3 F o g Mj8 ] W + aD l [ d4 F o c a \y} p }(x)=0 (3) 

for all x G [a,b], where F = L + Xg. 

Proof. Given hi,h 2 G V, |£i | <C 1 and \£ 2 \ <C 1, consider 

rb 

j(ei,e 2 )= / (Lo c a \y + e l h 1 + Bih 2 ]p){x)dx 
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rb 

z(ei,e 2 )= / (go%\y + e\hi + e2hz]p){x)dx-l. 

J a 

Since y is not an extremal for /, there exists a function ho such that 




and by the implicit function theorem, there exists a C 1 function £2(-)> defined in 
some neighborhood of zero, such that 

z(£l,£ 2 (£i))=0. 

Applying the Lagrange multiplier rule (see, e.g., [33, Theorem 4.1.1]) there exists a 
constant X such that 

V(;(0,0)+Ai(0,0)) = 0. 
Differentiating j and ; at (0,0), and integrating by parts, we prove the theorem. 

Example 1. Let y(x) = E a (x a ), x e [0, 1], where E a is the Mittag-Leffler function. 
Then q D" [y] = y. Consider the following fractional variational problem: 

J(y) = £(CD?\y}(x)) 2 dx^extx, 

i(y)= [ 1 y(x)$D?\y}(x)dx = l, 
Jo 

y(0) = l, y(l)=y u 

with / := / (y(x)) 2 dx and y\ := E a (l). In this case function F of Theorem 4 is 

Jo 

F(x,y, v, w) = v 2 + Xy{x) v 

and the fractional Euler-Lagrange equation (3) is 

x D«[2 C My]+Xy](x)=0. 

A solution to this problem is A = —2 and y(x) = y(x), x € [0,1]. 

The case when y is an extremal of / is also included in the results of [8]. 

Theorem 5. Ify is a minimizer or a maximizer of J on S, subject to the isoperimetric 
constraint I(y) = I, then there exist two constants Xq and X, not both zero, such that 

(**°SM/0 W+xD? [3*°SM/»] (x)+J)P [d 4 Ko c a \y] p ] (x)=0 
for all x G [a,b], where K = XqL + Xg. 
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Proof. The same as the proof of Theorem 4, but now using the abnormal Lagrange 
multiplier rule (see, e.g., [33, Theorem 4. 1.3]). 



5 Transversality conditions 

We now give the natural boundary conditions (also known as transversality condi- 
tions) for problems with the terminal point of integration free as well as yi,. 
Let 

& := {(y,x) e C l ([a,b}) x [a,b] \y{a) =y a ) . 
The type of functional we consider now is 

J(y,T)= [ T (Lo c a [y])( x )dx, (4) 

J a 

where the operator ^ [•] is defined by 

SM:=*^M*):a^M(*))- 

These problems are investigated in [1] and more general cases in [2]. 

Theorem 6. Suppose that iy,T) £ & minimizes or maximizes J defined by (4) on 
Then 

(d 2 Lo C a [y}) (x)+ x D$ [d 3 Lo C a [y]] (x) = (5) 
for all x £ [a,T]. Moreover, the following transversality conditions hold: 

(LoC M )(r)=0, x I 1 T - a [d 3 Lo c a [y]](T)=0. 

Proof. The result is obtained by considering variations y + eh of function y and 
variations T + eAT of T as well, and then applying the Fermat theorem, integration 
by parts, Leibniz's rule, and using the arbitrariness of h and AT. 

Transversality conditions for several other situations can be easily obtained. 
Some important examples are: 

• If T is fixed but y{T) is free, then besides the Euler-Lagrange equation (5) one 
obtains the transversality condition 

J l T - a [d 3 Lo c a \y]] (T) = 0. 

• If y(T) is given but T is free, then the transversality condition is 

(LogM) (T) -y\T)-J l T - a [d,Lo c a [y]] (T) = 0. 

• If y(T) is not given but is restricted to take values on a certain given curve i.e., 
y(T) = yf(T), then 
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(V(T)-y'(T)) -J l T - a [d 3 LoC[y]] (r)+ ( Lo c M) (r) = 
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